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Abstract 

In this paper we present a class of higher derivative theories of gravity which admit Birkhoff's 
theorem. In particular, we explicitly show that in this class of theories, although gcnerically the 
field equations arc of fourth order, under spherical (plane or hyperbolic) symmetry, all the field 
equations reduce to second order and have exactly the same or similar structure to those of Lovelock 
theories, depending on the spacetime dimensions and the order of the Lagrangian. 
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1 Introduction 



In General Relativity, Birkhoff 's theorem states that the spherically symmetric solutions of Einstein's 
equations in the vacuum are locally isometric to the Schwarzschild solution. This theorem has recently 
been generalized for a class of higher curvature theories namely the Lovelock theories [T] . These theories 
are natural generalizations of Einstein's theory in higher dimensions. While the field equations of a 
general higher curvature theory involves second derivatives of the Riemann tensor, Lovelock theories 
share the property of Einstein gravity that no derivative of the curvature tensor arise and hence the 
field equations are second order in the derivative of the metric [2J. In D-spacetime dimensions the 
general Lovelock Lagrangian consists of an arbitrary linear combination of all the 2fc-dimensional 
Euler densities, where 2k < D. Birkhoff's theorem in Lovelock gravity states that spherically (plane 
or hyperbolic) symmetric solutions of the Lovelock field equations are isometric to the corresponding 
static Lovelock black hole solution. The proof relies on the second order nature of the Lovelock field 
equations. The admittance of Birkhoff's theorem is related to the lack of spin-0 mode excitations 
in the linearized field equations. It is of fundamental importance in proving uniqueness theorems in 
General Relativity. 

Lately, there has been a great interest in other higher curvature theories of gravity coming from 
high energy physics, cosmology and astrophysics. These theories generically have higher order field 
equations in the metric formalism. However, as of yet, none of these theories are known to admit 
Birkhoff's theorem. In this article, we present a class of higher derivative theories which includes the 
Lovelock theories as a subclass and admits Birkhoff's theorem. These theories are characterized by the 
Lagrangian densities which give second order traced field equations and further all the field equations 
reduce to second order for spherically (plane or hyperbolic) symmetric ansatz. 

2 Construction of Lagrangian densities 

It is well known [3] that in dimensions greater or equal to four, the most general Lagrangian density 
which is quadratic in the curvature and gives second order traced field equations can be expressed as 
an arbitrary linear combination of the Gauss-Bonnet density and the quadratic conformal density 0. 
The Gauss-Bonnet density being the quadratic Lovelock Lagrangian gives second order field equations 
for generic metrics. The field equations obtained from the quadratic conformal density on the other 
hand is generically of fourth order. However, the trace of the field equations is proportional to the 
invariant itself, the proportionality factor being (4 — D)/2. This is not very surprising if one considers 
the variation of the action under infinitesimal conformal rescalings of the metric 5g ab = Qg ab which 
gives 

81 = J V=g£ ab 5g ab = J 5{^-gC) = J Q (2 - ^) y^A (1) 

1 Hereafter, by a conformal density we mean a scalar invariant constructed by contracting all the indices of a certain 
number of conformal tensors among themselves times yj—g. 
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where C = C ab C cd . It turns out that in dimensions greater or equal to four the Gauss-Bonnet 
density and the quadratic conformal density are the only two linearly independent invariants which 
generically gives second order traced field equations. Thus there is a two dimensional space of quadratic 
invariants with this special property in dimensions greater or equal to four. 

Now let us consider the same question in three dimensions. One might naively be tempted to think 
that since in three dimensions, the conformal tensor identically vanishes and so does the Gauss-Bonnet 
invariant, there are no quadratic invariants which gives second order traced field equations. However, 
this assertion can be explicitly shown to be incorrect. One may easily check that the density 

(2) 

does give field equations whose trace is of second order and is proportional to the density itself 0. In 
fact this is the only quadratic invariant which has this special property in three dimensions. Let us now 
ask ourselves the same question for cubic densities. First let us try to guess all the linearly independent 
cubic invariants which has this property. Obviously, first there is the six-dimensional Euler density all 
of whose field equations are of second order. Next one can repeat the same argument for the quadratic 
conformal density to the cubic conformal densities. Note that in dimensions greater than five there are 
two independent ways of contracting three conformal tensors. Hence, in dimensions greater or equal 
to five, there are at least three independent invariants which give second order traced field equations, 
namely- the six-dimensional Euler density and the two independent cubic conformal densities. In [5], 
we had explicitly shown that these are the only invariants which possess this property. However, in 
lower dimensions the Euler density identically vanishes and the two conformal densities are not linearly 
independent. In fact, in four dimensions the only cubic invariant which has the property is given by 
the conformal density. Whereas in five dimensions in addition to the cubic conformal density there is 
another independent density which shares the property. It is given by 

9/1 7? °d r? be r> a , r> cdr> ab r> , /in J? cd n c n d , OZl ' J n b r> cr> a z?&r?oi?_i_ r?3 

MK ab K cd K e + —K ab K cd K + WK ab K a K b + —K a K b K c - —K a K b K+—K . 

(3) 

This invariant is the cubic counterpart of the special quadratic density ([2]) in three dimensions. Now, 
to generalize the case for arbitrary higher order one needs to classify or understand the nature of these 
special invariants. 

Let us scrutinize the special quadratic invariant more carefully. First realize that in dimensions four 
and higher there are three linearly independent Riemannian invariants namely - R abcd R a b c d, R ab R a t, 
and R 2 . However, in three dimensions only two of them are independent. In other words any one 
of the three invariants can be expressed in terms of the other two. This identity is responsible for 

2 This density supplemented by an Einstein-Hilbert term and a cosmological constant has lately gained a lot of 
attention in the community as a toy model for quantum gravity in three dimensions and is known as the theory of New 
Massive Gravity [4] 
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vanishing of both the Gauss Bonnet density and the conformal density Analogously, there are eight 
linearly independent cubic Riemannian invariants in dimensions greater or equal to six. However, in 
five (and lower) dimensions the six-dimensional Euler density vanishes identically. In fact, this is the 
only independent identity in five dimensions among the cubic scalar invariants. So, this identity is also 
responsible for the linear dependence of the two conformal invariants. The same occurs for invariants 
of arbitrary higher order k in dimensions D = 2k — 1 because of the following relation 

£C\d\---c k d k r> aifei . . . E> a k b k _ rcidi— cj,d k s~< a\bi p a k b k _ q /a\ 
a^bi-a k b k cidi c k d k a\bi---a k b k cidi c k d k ' 



where 6..'. is the generalized Kronecker delta. Now, in the case of quadratic invariants, when one takes 
a particular linear combination of the Gauss-Bonnet density and the conformal density in arbitrary 
dimensions and reexpresses the conformal invariant in terms of the Riemannian invariants then it 
factorizes by (D — 3). Explicitly, 

R abcd Rabcd ~ *R ab Rab + R 2 ~ C abcd C abcd = [AR ab R ab - ^y# 2 ) ■ (5) 

Obviously, since the left hand side is a linear combination of the Gauss-Bonnet density and the 
conformal density, it gives second order traced field equations in arbitrary dimensions. This implies 
that the invariant inside the parenthesis on the right-hand side also gives second order traced field 
equation in arbitrary dimensions. However, the left-hand side vanishes in three dimensions whereas 
the term inside the parenthesis on the right-hand side does not. In fact in three dimensions it gives the 
special invariant mentioned earlier. Similar situation arises in the cubic case where a particular linear 
combination of the six-dimensional Euler density and the two conformal invariants can be factorized 
by D — 5 and the remaining invariant does not vanish in five dimensions identically. This gives the 
special cubic invariant in five dimensions ([3]). Interestingly this can be generalized to arbitrary higher 
orders when one realizes that there is always a particular combination of the 2fc-dimensional Euler 
density and the fc-th order conformal invariants which can be expressed as 

?cidi---c k d k / p aifei p a k b k _ ^ a-^bx (~i a k b k \ 

°a 1 b 1 -a k b k \ rL c 1 d 1 '" a c k d k U cidi "'Vt >> W 

which when expanded (reexpressed) in terms of Riemannian invariants is factorized by (D — 2k + 
1). The remaining invariant is non- vanishing in dimensions D = 2k — 1 but vanishes identically in 
lower dimensions. Based on this observation we had proposed a conjecture in [5] that any kth order 
Riemannian invariant which gives second order traced field equations can be expressed as y/—g times 
a linear combination of the special invariant, all the linearly independent conformal invariants and a 
divergence term@. Explicitly, if in D dimensions, there are linearly independent k(> 2)th order 
conformal invariants for {i = 1, • • • N^}, then the most general action which gives second order 



3 Of course this can be equivalently understood in terms of vanishing of the conformal tensor. However, here we are 
concerned with scalar identities. 

4 We have checked this conjecture up to order four. 
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traced field equations is given by 



where 



( 



N 



(fc) 
D 



(Xq M^ 1 + ^ + a divergence term 

D - 2 



i=l 

rci(ii---c fc (i fc / D 0161 



; / tj aiOi d a fc 6 fe 



c 



cidi 



Cfcdfc 



(7) 



(8) 



Note that in dimensions D > 2k, the special invariant J\f^ can be expressed as a linear combination 
of the 2/c-dimensional Euler density and all the conformal invariants and in dimensions D < 2k — 1, it 
vanishes identically. 

Now we show that a subclass of these theories which generically gives fourth order field equations 
but for spherically (plane or hyperbolic) symmetric spacetimes all the field equations reduce to second 
order!. And it is in this class where Birkhoff 's theorem can be shown to hold. 



3 Field equations for spherically (plane or hyperbolic) symmetric 
spacetimes 

Consider the general spherically (plane or hyperbolic) symmetric spacetimes given by the following 
line element 

ds 2 = g ij (x)dx i dx j + e 2X{x) dT, 2 , (9) 

where dY> 2 = g a p{y)dy a dyP is the line element of a (D — 2)-dimensional space of constant curvature 
7. Let V be the Levi-Civita connection on the two-dimensional space orthogonal to the constant 
curvature space and R be the corresponding scalar curvature. Then the nontrivial components of the 
Riemann curvature tensor and the conformal tensor are given by 

v fc = lm r 4 x = = ( io ) 

p ik \P_ ~ 3)3 rik p flX & irfiX p ifJ- ~ H XfM (llN 

jl "2(L»-1) J " vp " (D - 1)(D - 2) i" 2{D- l)(D-2) 3 v ' 1 ' 



where 



fi = 7 e- 2A -(V m A)(V m A), (12) 
A) = V*VjA + (V\\)(Vj-A), (13) 
and S = R + 2V fe V fe A + 2 1 e' 2X . (14) 



Note that since all the components of the conformal tensor are a mere multiple of the function S, 

(k) j—i ~h (k) 

each of the conformal densities W^> s evaluated on the metic (J9J) are proportional to S . Let Wm = 

uj m (D,k)S k . Then the field equations for the action ([7j) evaluated on the metric ansatz ([9]) are given 
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by 



( k ) i5fc-i 



, {k)i = (D - 2)l(D - 2)a), K> B 
' J 2{D - 2k - !)!(£> - 2k + 1) L 



2k5f l A{ — (D — 2k — l)8)B) 



N 



(k) 



+ k | a$u m (D,k) 

m=l 



(D-2)a {k) u (D,k) 
2 k (D-2k + l) 



v l J (s k - 1 ) 



(15) 



(D - 2)\af ] B k - 2 



y P~ (D-2k-l)\(D-2k + l)P 



(D - 2k - l)(D - 2k - 2)B 2 



+ k(R - 2(D -2k- l)Ai)B + 2k(k - l)ffji%A k 



( n d ] 

^m=l 



(fe ),., (D-2)a ( k) u (D,k) 



2 k {D-2k + l) 



gW = gW? = o 



5 a Q{S 



fc-l^ 



(16) 
(17) 



where V l j and Q are two (related) second order linear differential operators denned on the two- 
dimensional space orthogonal to the constant curvature base manifold and uJo(D,k) is a positive 
number (see Appendix). Notice that the fourth derivative terms arise when the operators Vj and Q 
act on the function S k ~ l . However, in all the field equations these terms are multiplied by a numerical 



factor which depends on the coupling constants a[ , for {i = 0, • • • N^'}, the dimensions D and the 
order k. Hence, if one chooses the coupling constants in such a way that this factor vanishes then the 
field equations reduce to second order. In particular, choosing {a^ ^ 0, a^} such that 



N 



(k) 



*£a$u; m (D,k) 



m=l 



(D-2)a { k) oo (D,k) 
2 k (D-2k + l) 



(18) 



give the same equations as in pure Lovelock gravity theories in dimensions D > 2k. Note that for 
any dimension D > 2k for a given k, = n!$ whereas n!£) = ^vjft—i + 1 as explained previously 
due to the identity (|J|). Therefore, the number of independent densities of order k > 2 satisfying (p~8|) 
in dimensions D > 2k is N^} ■ Also, note that in D = 2k or D < 2k — 1, all the field equations 
corresponding to the densities satisfying (|18p vanish identically when evaluated on the metric ([9|). 



4 Classification of theories admitting Birkhoff 's theorem 

Let us now discuss the different cases up to the first few orders. Since there are no conformal invariants 
of order one, the only non-trivial action, in any dimension D > 2, is given by the Einstein-Hilbert 
term. Next, there is only one conformal invariant of order k = 2 in dimensions D > 4. This indicates 
that there is only a one-parameter family of densities in D > 4 which admit Birkhoff's theorem. 
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So, this must be the Gauss-Bonnet density. One can explicitly check that this is indeed the case by 
calculating uiq(D,2) and u)\(D,2). In D = 4, the Gauss-Bonnet density is a topological term and in 
D < 4 it vanishes identically and hence there are no densities quadratic in curvature which admits 
Birkhoff's theorem in dimensions D < 40. Now, let us consider the cubic invariants. It turns out that 
there are two independent conformal invariants in dimensions D > 6. Therefore, there must be two 
linearly independent densities which admit Birkhoff's theorem in D > 6. Obviously, a particular linear 
combination of these two densities gives the third order Euler density, which in turn has second order 
field equations for any metric and is already known to admit Birkhoff's theorem. This implies that 
any other linearly independent combination of the two densities represent a four-derivative theory, 
whose field equations when evaluated on the metric (0) are either the same as that of cubic Lovelock 
theory or are identically vanishing. Whereas in D = 6, since the third order Euler density is a 
topological term, there is only one linearly independent non-trivial density which represents a four- 
derivative theory. However, when evaluated on the metric Q, all the corresponding field equations 
vanish identically. Now, in five dimensions, as explained previously, there is only one independent 
cubic conformal invariant. Also, the third order Euler density vanishes identically in D = 5. So, there 
is a unique cubic theory in D = 5 which is a four-derivative theory and admits Birkhoff's theorem. 
Even though there is no cubic Lovelock theory in five dimensions, the field equations of this theory 
when evaluated on the metric has a similar structure as that of the cubic Lovelock theory in higher 
dimensions. This theory was first presented in [7] and the Birkhoff's theorem was proven explicitly [§. 

So for arbitrary order k > 2, in dimensions 

(k) 

D > 2k: there are N^J independent densities whose field equations are generically of fourth order 
but when evaluated on the metric ([9]) either give the same field equations as those of the corresponding 
A;th order pure Lovelock theory or vanish identically. 

(k) 

D = 2k: there are AQ> independent densities out of which one particular linear combination gives 
the 2fc-dimensional Euler density which is a topological term. The field equations from any other 
linearly independent combination are generically of fourth order but when evaluated on the metric (JH) 
vanish identically. 

(k) 

D < 2k — 1: there are — 1 independent densities whose field equations are generically of 
fourth order but when evaluated on the metric Q vanish identically. This is because both the special 
invariant (JHJ) and the right hand side of the equation (|18p vanish. 

(k) (k) 

D = 2k — 1: there are = N^ k — 1 independent densities whose field equations are generically 

of fourth order but when evaluated on the metric Q either vanish identically or reduce to second 
order and has a similar structural form as that of pure kth. order Lovelock theory. In the later case 
the Birkhoff's theorem can be proved in the following way. 

Let us begin by fixing coordinates in the two-dimensional spacetime with metric gij, and the gauge 

Interestingly, a weaker version of Birkhoff's theorem has been shown to hold for conformal gravity in four dimensions 
[6|. There it has been proved that the most general spherically symmetric solution of conformal gravity is static up to a 
conformal gauge transformation. 

6 It was simultaneously found in [8] where it was named Quasi-topological gravity. 
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freedom by choosing Ve A ( x ) to be a spacelika^ vector on gij\ 



The equation Q t = implies 



while Q r = reduces to 



9ij- 

ds 2 = - fit, r)dt 2 + — r + r 2 d£ 7 . (19) 

9{t,r) 



( 7 - gf' 1 (2 ( 7 - g) + krg') = , (20) 



( 1 -g) k ~ 1 (2( 7 -g)f + krgf')=0, (21) 



where the prime on g and / denote a partial differentiation with respect to r. Focussing first in the 
non-degenerate case, i.e. g (t, r) ^ 7, from equation ([2UJ) we obtain 

g(t,r)=F 1 (t)r^ k + 7 , (22) 

where F\ (t) is an arbitrary function of t. Replacing this in equation (|2ip . we obtain 

f(t,r) = F 2 (t)g(t,r) . (23) 

The arbitrary function F2 (t) can be reabsorbed by a coordinate transformation without any loss of 
generality. Consequently 

f{t 1 r)=g{t,r)=F l {t)r 2 l k + 1 . (24) 
Equation = then implies dF\/dt = 0, which in turn implies that F\ is a constant c. Then 

(k)a 

the equations along the base manifold Q a = are satisfied without any further restriction, and we 
obtain the following metric 

ds 2 = _ ( cr 2/k + 7 ) dt 2 + -J^— + ^dT 2 , (25) 

where c is an integration constant. The spacetime is asymptotically locally flat. In the hyperbolic 
case (7 = — 1) the metric describes a black hole, provided c is positive. This completes the proof of 
the Birkhoff's Theorem for the non-degenerate case. 

In the degenerate case where g (t, r) = 7, all the equations are trivially satisfied, and the metric 
function / (t,r), is left undetermined. Birkhoff's theorem does not hold in this case. 

One can also consider actions which are non-homogenous in the order k. In such a case, unless there 
is a non- vanishing contribution to the action from the invariant ([8|) of order k = -^p- , the spherically 
(plane or hyperbolic) symmetric solutions are given by the corresponding solutions of general Lovelock 
theory of order k = [^F^] • However, if there is such a contribution (necessarily in odd dimensions) 
then the general spherically (plane or hyperbolic) symmetric solution is again static but does not 



7 The proof of the Birkhoff's theorem for timelike and null cases follows along the same lines. 
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belong to the family of solutions of general Lovelock theory. Nevertheless, the general solution has 
the same structural form as that in general Lovelock theories of order k = and can be implicitly 
written in terms of an algebraic equation for the metric function gu- Generically, the solution is of 
the form 



ds 2 



-f(r)dt 2 + ^- + r^ 
f(r) 



(26) 



where the function / (r), solves the following polynomial equation 



m 



£ (D - 2k) r 



D~2k-1 



(7 -/(or 



M 



(27) 



fc=0 



where M is an integration constant and 



(fc) _ (D-2) (D-2)\ (fc) 
°° ~ (£>-2fc + l)! 



(28) 



Therefore, in any dimension D, in addition to the cosmological constant, there is a p-parameter 
family of (non-trivial) Lagrangian densities which generically admit Birkhoff 's theorem, where 



p=l + 



E * J 



+ 



\i<fe<-^ 



\ 



for even 



1 + 



E 



A/ 



(fc) 

2fc 



1 + 



Vi<fc<^ 



E w 



for odd D. 



Here, for both even and odd D, the unity on the right hand side corresponds to the Einstein-Hilbert 
term, whereas the second and the third term correspond to the terms in the action which generically 
gives fourth order field equations but when evaluated on the metric ansatz Q respectively reduces to 
second order and vanish entirely. 



5 Conclusions 

Our findings raise several natural questions. Firstly, is there an even wider class of theories which 
admit Birkhoff 's theorem? Secondly, we have seen that the solutions of the theories discussed here all 
have the same structural form namely those of Lovelock gravity theories. Is there any theory which 
admits Birkhoff 's theorem but the corresponding solution does not have the same structure? It seems 
to us that, if it is necessary for Birkhoff 's theorem to hold, that the field equations reduce to second 
order under spherical (plane or hyperbolic) symmetry then the corresponding solutions will have the 
same structure. 
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Appendix 

The operators and Q in (|15p and ()16|) are second order linear differential operators defined on the 
two-dimensional space orthogonal to the constant curvature base manifold and are given by 



' j 



S) [ I + [D - l)V fc V fc A + {D-2){D- l)V fc AV fc A + V k V k + (2D - 3)V k XV k - ^ 

(D - 2)(V i V j X + DV'XVjX) - V*Vj — (D — l)(V*AVj + VjAV') 

D 



a 1 



V % - - S 1 

2k 



D-2 

The numerical factor Uq(D, k) is given by 



(29) 
(30) 



MD, k) = — _ _ [(D - 2k)(D - 2k - 1) + k(k - 2)D{D - 3) + k(k + 1)(D- 3)] 

(31) 
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